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The discrepancy between the PQCD calculation and the CLEO data for xa — >■ jV (1^ = p", uj, 4>) 
stimulates our interest in exploring other mechanisms of Xci decay. In this work, we apply an im- 
portant non-perturbative QCD effect, i.e., hadronic loop mechanism, to study Xci — ^ ^V radiative 
decay. Our numerical result shows that the theoretical results including the hadronic loop contri- 
bution and the PQCD calculation of Xa — >■ yV are consistent with the corresponding CLEO data 

(^ • of Xci — > yV. We expect further experimental measurement of xa — i> yV, which will be helpful to 

C^ ' test the hadronic loop effect on Xci decay. 
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' _ In the past three decades, a series of the observations of S-wave, P-wave and D-wave charmonia make charnionium 

iT^. family abundant. Nowadays, charm physics is still an intriguing research field with challenges and opportunities [l|. 
d ' Especially, the study of charmonium may provide valuable information on non-perturbative QCD effects. 

As an important and effective approach to deeply learn the underlying properties of charmonium, charmonium 
decay is an extensively focused research topic. Among the observed charmonium states, J/ip and iIj{2S) are of 
abundant experimental information of decay just listed in Particle Data Group (PDG) [3|- However, the experimental 
vQ ' measurement relevant to the decay of P-wave charmonium is far less than that of J ftp and tjj{2S). Thus, more 
\^ [ experimental and theoretical explorations of P-wave charmonium decay are becoming active, especially with the 
O ' running of BES-III. 

^D ' Recently, the BES-III Collaboration announced its observations of Xcj{J = 0, 1, 2) decaying into two light vector 

mesons (3[. Among these decay modes of XcJ to two light vector mesons, the OZI suppressed processes Xci -^ wo;, cfxj) 
f^ ' and the double-OZI suppressed process Xci ~^ w0 were firstly observed. In order to explain the evasion of the helicity 
^^ . selection rule in these processes, the hadronic loop effect, an important non-perturbative effect relevant to the decay 
of charmonia [4j-|8| and molecular system pj|-(2l|. is introduced in Refs. J22h24| . which also indicate that applying the 
hadronic l oop mechanism to other Xcj decays will be helpful to further test the hadronic loop effect on Xcj decay. 

In Ref. [2^ , the radiative decays of charmonia J/ip and Xcj into light meson are studied by the perturbative QCD 
(PQCD) approach, where a complete numerical calculation for the quark-gluon loop diagrams was performed. The 
^ . obtained theoretical results for J/tp — ^ 777, 777' can well reproduce the experimental data. Furthermore, the branching 
ratios of Xcj -^ 7P°, 7^, l<t> were predicted, which arc B{xci -^ 7P°) = 1-4 x 10~^, B{xci —> 7^) = 1.6 x 10~^ and 
'S(Xci -> 70) = 3.6 X 10-6. 

In 2008, the radiative decays of Xcj were first measured by the CLEO Collaboration using a total of 2.74 x 10^ decays 
of the -0(25') collected with the CLEO-c detector [H]. The reported results are ;B(Xci -^ 7P°) = 243±19±22x 10"^, 
I3{xci -^ 7^) = 85 ± 15 ± 12 X 10"^ and B{xci -^ 7'/') = 12.8 ± 7.6 ± 1.5 x 10"^. One notices that the experimental 
results of the branching ratios of Xci ~^ IP^i 7^, 70 are an order magnitude larger than the corresponding theoretical 
predictions |25| . 

The above difference between the experimental measurement from CLEO and the theoretical result calculated in 
PQCD shows that there should exist the extra effect on the Xci radiative decays into a light vector meson, which 
stimulates our interest in exploring the underlying mechanism to resolve this large discrepancy between the theoretical 
prediction by PQCD [2^ and the experimental measurement by CLEO [20] . 
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the hadronic loop effect can explain the experimental observation of the OZI 
p and the double-OZI suppressed process Xci -^ ^4* well. In this work, we extend 



As indicated in Refs. [2J, [24 
suppressed processes Xci — > wo;, (jjcj) 
the hadronic loop effect to xa -^ jV {V = p^, oj, cj)) process to answer whether large discrepancy between the 
theoretical calculation and the experimental measurement of Xci -^ "fV can be alleviated. Xci ~^ 7^ is similar to 
Xci -^ yy since both processes occur via the intermediate charmed mesons under considering the hadronic loop 
effect. What is more important of this work is to offer an effective approach to test the proposed hadronic loop effect 
applied to explain Xci -^ yy decay processes. 

The paper is organized as following. After the introduction, we present the formula of the hadronic loop contribu- 
tions to Xci radiative decays to a light vector meson, which includes the effective Lagrangian employed in this work 
and the decay amplitudes. In Sec. IIIIl the numerical results of xa -^ IP^, 7^, 70 are given. The last section is the 
discussion and conclusion. 



II. HADRONIC LOOP EFFECT ON xa -> -yV 



As indicated in Ref. (23|, the nonperturbative QCD mechanism, i.e., hadronic loop effect, may play a crucial role in 
understanding Xci decay. In Table. [TJ the typical diagram depicting the hadronic loop effect on Xci -^ 7^ at quark 
level is given, which is different from the quark-gluon loop diagrams depicting x^i ~^ 7^ at PQCD approach in Ref. 
[25| . The transition element of Xci — >■ 7^ can be expressed as 



Mixci ^ ly] = J2(^y\n^'^\i){i\H^'^\xci), 



(1) 



which reflects the intermediate state contribution to Xci -^ 7^- Here, H^^' represents the interaction of Xci and 
DD* + h.c. and TiS^i describes interaction DD* + h.c. — > ^V by exchanging an appropriate charmed meson. 



Quark level 



Hadron level 




FIG. 1: The quark level typical diagrams (the first column) describing the hadronic loop efi'ect on Xci — > 'yV and the hadron 
level schematic diagram corresponding to Xci — >■ 7p'^ (the second column). The red and black lines denote the charm quark and 
light quark. The photon emits from charm quark line or light quark line. By the charge conjugate transformation !)(*)+ ^ £)(*)o 
and D^*' ^ D^*'^, the rest two diagrams of Xci — >■ DD* + h.c. — )> ^p^ can be obtained by diagrams (a) and (d). 



With Xci — >■ 7/3° as an example, we list the schematic diagrams at hadron level in Table [TJ where Xci first dissolves 
into two virtual charmed mesons which is originated from the coupled channel effect. Then these two virtual charmed 
mesons DD* + c.c. turn into a photon and p" meson by exchanging the charmed meson. Due to the mass of Xci being 
lower than the threshold of DD* , the charmed mesons in the loop are off-shell. 

In the following, we still use Xci ~^ IP^ as an example to illustrate the relevant calculation of hadronic loop diagrams 
listed in Fig. [U where the effective Lagrangian approach is applied to write out the decay amplitude throughout this 
work. 

The Lagrangian for Xci2?2?* coupling reads [23] 



^xi-DV = igxaiV-D'Xci ■ 2?i 2?' + h.c. 



(2) 



The effective Lagrangians responsible for 2?(*)p(*)y interactions are 

where !?(*) = (£>(*)",£>(*)+, D,l*'+) and A d ^, B = A{d^,B) - {d^A)B. The matrix of the nonet vector mesons V is 
defined as 



V = 






(4) 



We need to specify that the Lagrangians in Eq. (j3|) is just the first term in an infinite series of terms that represents 
the hadronic representation of the QCD Lagrangian, which is restored from the Lagrangian constructed in the chiral 
and heavy quark hmits J28l - l32| . Thus, we adopt usual z/(fc^ — rn^) and i{—g^^ + k^k'' /'m?)/{k'^ — m?) propagators 
when writing out the decay amplitude of Xci -^ DD* + h.c. -^ 7p°. The relevant coupling constant will presented in 
the next section. 

The Lagrangians for ^DD and ^D* D* interactions can be obtained from the Lagrangian for free scalar and massive 
vector fields by the minimal substitution d'^ ^^ d^ + ieA'^ , which are p33j 

C-D-D-^=ieA^,D- d D+ + ieAf,D^ d £>+, (5) 

Cw'j = ieA^ |.9"^^r a'' D*+ + g^t" D^' d'' D*+ - g'^«a^Z?r i?;^ 

MeA, |g"^i?:- T D*+ + g^PD^-d^D:; g'^-d^D^-Di; \ , (6) 

respectively. Here, the electromagnetic interactions of D^D'^^ and D*'^D*'^'y does not exist. 
The Lagrangian describing the electromagnetic D*Dj vertex is [Sj] 

Cv'V^ = S^lgD^^D^.e^'^^F^^Dl+D- + ^go^o^o^e^'^-PF^^VllD^X + h.c . (7) 

where Ff^^, — d^A^ — d^A^ and D*^j^ = 9^iZ?*"'+ — 9^1?*°'+ is the stress tensor of the vector charmed meson. To 
some extent, the Lorentz structure just shown in Eq. ^ is same as that in Ref. [3J, |35|- In Eq. ((T]), parameters 
5d*+d+7 and gD*°D°f a-rc introduced to get consistent results with experimental measurements of D*+ — ^ Z)+7 and 
£)*o ^ £)0^_ rpi^g theoretical decay widths of D*+ -^ D+-f and D*^ -^ D^j are 

rp*+ ^ Z?+7) = Y^9Un.,<^. (l - ^) , (8) 

r(D*o ^ D'j) ^ ^^gl,.^o,ml,. (l ~ ^) . (9) 

According to the experimental widths r{D*+ -^ D+j) = 1.54 keV [1] and r(L>*" -^ D"j) = 26.04 keV [2, [13, the 
coupling constant go'D-y is fixed as 

|5^.+^+^|=0.5GeV-\ |to-"DO-,| = 2.0GeV-^ (10) 

Both calculations based on Lattice QCD [401 and QCD sum rules (QSR) [35| predict that the coupling constant for 
the radiative decay of the neutral charmed meson has a positive sign while the coupling constant for the charged 
charmed meson radiative decay is negative. In present work, we follow such a convention and take go'+D+'y = —0.5 
GeV~^ and gD''>D«'y = 2.0 GeV~^. For the coupling constant of D*Dsj interaction, the calculation from QSR gives 
9D*Ds-y = —0.3 ± 0.1 GeV~^ [33 ■ In present work, the central value is adopted. 



According to the Lagrangian just listed above, we obtain the decay amphtudes of Xci ^^ IP^ corresponding to the 
diagrams in Fig. [TJ 
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which are resuhed from the charge intermediate charmed mesons. Thus, the total decay amplitude of Xci -^ DD* 
h.c. — > 7p° is 
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(15) 



where subscripts C and N denote the corresponding amplitudes being from charge charmed meson loop and neutral 
charmed meson loop, respectively. Ai^ and Ai}^ is obtained by amplitudes ^A^ and Ai^ with the replacements 
of the mass and coupling constants, i.e., gD'+D+-y ~^ gD*°D°li "?,£)(.)+ -^ rrij^Mo and mjj{,)- -^ 77i£i(.)o. In Eqs. 
(fTT|) -(|14 p . the form factor J^{q^) is introduced to depict the inner structure of the interaction vertex of the exchanged 
charmed meson and the intermediated state. As what we have done in Ref. [23|, a dipole form of the form factor is 
employed 



Hq') = 



A2 



A2-g2 



(16) 



Furthermore, the form factors with the pole form also play the role to make the ultraviolet divergence disappear, 
in analog to the cutoffs in the Pauli- Villas renormalization scheme. Here, the cutoff A can be parameterized as 
A = TO + aAqcD with Aqcd = 220 MeV and to is the mass of the exchanged meson in Fig. [T] [S^ ■ We emphasize 
that a dipole form factor is introduced in the numerical calculation of this work, which was applied to the calculation 
of Xci -^VV {V = p, w, (j)) in Ref. JH. 

The decay amplitude of Xci — > DD*+h.c. -^ 70; is similar to that of Xci — ^ DD*+h.c. — !> -fp^ besides multiplying an 
extra factor -1 since the DD*V vertex would have a different sign for DD*ui and DD* p due to 5't/(3) flavor symmetry. 
In addition, we need to replace the mass of p with those of uj. For Xci -^ DsD* + h.c. -^ jip, the corresponding 

charmed -strange mesons. The total decay amplitude is 



hadronic loops are composed of Dg* and Dg* 
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Mixci -^ DsD* + h.c. -^ 70) = [M'^c' + ^c' + ^c' + -^c 



(17) 
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with the replacements of the corresponding masses and coupling constants. 



III. NUMERICAL RESULT 

Before performing the numerical calculations, we need to introduce the coupling constant relevant to the effective 
Lagrangian listed in the previous section. The coupling constants for XciT^T^* smd VDV interactions are listed in 
Table, n 



Coupling constant 


Expression 


Value 


Coupling constant 


Expression 


Value 






-21.44 GeV 


SXalOsD* 






ffXclDD* 


■Z^-Zgi^mDmD'm^^-^ 


■z^ ■Agi^mD,mD*m^^^ 


-22.60 GeV 


g-D-DV 


Pgv/V2 


3.71 


g-D'-D'V 


Pgv/V2 


3.71 


f-D'-DV 


\gv/V2 


2.31 GeV"' 


fo'D'V 


\mD*gv/V2 


4.64 


fo^D'V 


\Tno%gv 1 ^ 


4.88 









TABLE I: The coupling constants relevant to the calculation of xa ^^ 'yV. Here, gi is related to the XcO decay constant f^^o 
|£S--p— with /x^o - 0.51 GeV [sj]. Other parameters include gv = mp/f^, nip = 0.77 MeV, /3 = 0.9, 



via relation gi = 

A = 0.56 GeV-\ g = 0.59 and /^ = 132 MeV [32, H H [39 



With the above preparation, the radiative decays of Xci to a light vector meson are estimated. As a free parameter, 
a is introduced in the cutoff A of the form factors, which is usually dependent on the particular process and taken to 
be of the order of unity. In Fig. [2l we present the branching ratios of Xci — ^ IP^i 7^, 70 dependent on the parameter 
a. For comparing with the experimental data [2g], the theoretical result includes the hadronic loop contribution 
obtained in this work and the PQCD estimation in Ref. [25| . 
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FIG. 2: (color online) The branching ratios of xa — >■ 7/9 , Ji-^, ')4> dependent on the parameter a. The red dashed-lines with 
the blue bands are the experimental measurement [M. The blue solid lines correspond to the theoretical calculations including 
the hadronic loop effect and the PQCD calculation [2^. The vertical yellow bands in the sub-figure denote the overlap of our 
results with the corresponding experimental measurement. With the same a range sandwiched between two green vertical solid 
lines, the obtained branching ratios of xa — ^ 7p'', 71^, 7<^ in this work are consistent with the experimental data. 



As shown in Fig. [2l there exists overlap between the numerical result obtained in this work and the experimental 
data announced by CLEO. The corresponding a ranges for Xci -^ 1P°, 7W, 71^ are 2.18 < a < 2.35, 2.06 < a < 2.28 
and 1.16 < a < 2.77, respectively, which are in the reasonable parameter space. Especially, we need to emphasize 
that there exists a common a range 2.18 < a < 2.28 for Xci -^ 7P°, 71^, 70 radiative decays. 



IV. DISCUSSION AND CONCLUSION 



As an important non-perturbative QCD effect, the hadronic loop mechanism was proposed in studying J/ip and 
■0(3770) decays [J, la, Q- By using this mechanism, the OZI suppressed processes Xci -^ W with VV = wo;, 00 and 
the double-OZI process Xci "> "^0 announced by BES-III were explained well in the recent work j23Ll24|. 

The CLEO Collaboration announced the experimental results of Xci ~> "fV m 2008 [2a], which are an order 
magnitude larger than the corresponding theoretical estimations calculated by PQCD [25[. To search for the source of 



the discrepancy between the PQCD calculation [25| and the CLEO data [26[ for Xci — > jV, in this work we propose 
the hadronic loop contribution to Xci — ^ jV- Under the hadronic loop mechanism, Xci — ^ iV is similar to Xci -^ VV 
process, both of which occur via the intermediate DD* just shown in Fig. [T] To some extent, the study of this work 
can be as a test for the hadronic loop effect, which was applied to explain Xci ~^ VV decays [23ll24|. 

Our numerical result of Xci ~^ iV indicates that the theoretical result including the hadronic loop contribution 
and the result in PQCD calculation for Xci -^ iV can reach up to the experimental data of Xci ~^ iV . Thus, non- 
perturbative QCD effect, i.e., hadronic loop mechanism, can be as the underlying source to alleviate the difference 
between the PQCD calculation and the CLEO data of Xci -^ iV . As indicated in Refs. I23,l2j], the hadronic loop 
effect also plays an important role to Xci -^ VV. Thus, the success of explaining xa — >■ VV and xa -^ iV under 
the hadronic loop mechanism not only tests the model itself, but also shows that the non-perturbative effect on Xci 
is important. Further experimental and theoretical studies of Xcj decay are encouraged. 
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